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ABSTRACT 



The density of the far infrared / submillimeter (FIR/SMM) diffuse 
extragalactic radiation field has recently been determined from COBE data. 
Nearly simultaneously, deep FIR/SMM surveys have detected substantial 
numbers of optically unidentified sources, which have led to the proposal 
that galaxies and protogalaxies at red shifts z = 2 to 4 may account for an 
appreciable fraction of the background. Here I show that, if the reported 
radiation levels are generated through nucleosynthesis, most of this energy must 
have been produced at epochs z ^ 2. Hubble Deep Field data cited by Madau 
et al. (1998) indicate that the bulk of the integrated extragalactic background 
must have been generated even more recently, at z < 1. 
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1. Introduction 

Two critical observational advances have taken place in the past two years. The first 
is the determination of a diffuse extragalactic FIR/SMM background (Puget et al. 1996; 
Hauser et al. 1998). The second is the detection of faint distant galaxies at FIR/SMM 
wavelengths (Kawara et al. 1998; Hughes, et al. 1998; Barger et al. 1998). Hughes et 
al., who found no visible or infrared galaxies at the positions of most of the sources they 
reported in the Hubble Deep Field, suggest that these and much of the SMM background 
originate in distant, dust-shrouded galaxies at red shifts z = 2 to 4, undetectable at short 
wavelengths. 

A number of observations argue against such an early origin for most of the diffuse 
background. Lanzetta et al. (1995) have suggested that most Ly-a absorption systems are 
progenitors of galaxies, hydrogen rich at high red shifts and rich in stars at low red shifts. 
This is supported by Rao, Turnshek, and Briggs (1995), who found a drop in gas content by 
an order of magnitude between z ~ 3.5 and 0, but roughly comparable density parameters, 
for gas at high red shifts, fl g (z = 3.5), and luminous matter (stars) seen today, fl s (z = 0). 

Pettini et al. (1998) have traced the evolving abundance of zinc in damped Lja systems 
from red shifts z = 3.4 to 0.4. Zinc is a useful tracer of metallicity, since it is not readily 
deposited from the gaseous phase onto interstellar grains. At red shifts beyond z ~ 2.8 
only upper limits could be established, indicating a zinc deficiency > 25. But even at red 
shifts z = 1 the zinc abundance remains a factor of ~ 3 to 10 below the solar abundance. 
At red shifts z ^, 3, the abundance of carbon in these systems is roughly a factor of 100 
lower than solar, and independent of the optical depth (Songaila & Cowie 1996; Cowie & 
Songaila 1998). These results do not appear to be seriously affected by dust absorption, 
since the extinction even for blue light reaches unit optical depth only rarely, at the highest 
column densities, ~ 2 to 3 x 10 21 cm~ 2 (Zuo et al. 1997, Songaila et al. 1995). The cited 
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abundances of heavy elements thus are likely to be reasonable tracers of energy generation 
in massive and supermassive stars at early epochs. These metallicities also agree with data 
for the oldest Galactic stars formed at about the same epoch (Timmes et al. 1995; Cowan 
et al. 1995). 

The general picture that emerges is that the earliest production of heavy elements 
occurred in Population III stars whose masses were sufficiently high to produce and 
explosively eject r-process elements. The metallicity at these epochs was 2 to 3 orders of 
magnitude lower than today. Substantial accumulations of heavy elements did not appear 
until z ~ 2. 

Sections 2 and 3, below, respectively summarize the background observations and 
FIR/SMM surveys. Section 4 discusses energy generation rates, while section 5 shows why 
the primary epoch of energy generation must lie in the relatively recent past, at red shifts 

2. Background Observations 

Using the Far-Infrared Absolute Spectrometer (FIRAS) onboard the Cosmic Microwave 
Explorer satellite (COBE), Guideroni, et al. (1997) estimated a flux at wavelengths 
between 200 and 400 /xm, of \l\ ~ 7nWm~ 2 sr^ 1 at ~ 300/im, corresponding to a 
radiation density p300Atm ~ 3 x 10" 15 erg cm" 3 . Measurements from the Diffuse Infrared 
Background Experiment (DIRBE) aboard the same satellite led Hauser et al. (1998) to 
a somewhat higher integrated energy from 140 to 240 jum, of \l\ ~ 10nWm _2 sr~ 1 , or 
P2oo M m ~ 47rA/ A /c ~ 4 x 10" 15 erg cm" 3 . Fixsen et al. (1998) used both DIRBE and FIRAS 
data and found the total FIR/SMM background from 125/jm to 2 mm to be slightly higher, 
Pfir/smm ~ 6 x 10~ 15 erg cm" 3 . 

These data represent only a fraction of the total diffuse extragalactic energy density. 
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In the mid-infrared, between ~ 5 and 100 //m, only coarse upper limits are available, all 
of which are much higher than the FIR/SMM values. At 3.5 /im, Dwek & Arendt (1998) 
obtain a radiation density of 3 x 10 -15 erg cm -3 from DIRBE observations. At 3600 to 
22,000 A, Pozzetti et al (1998) have estimated a background flux of ~ 4 x 10~ 15 erg cm" 3 , 
from discernable sources in the Hubble Deep Field. We can write the total extragalactic 
radiation density, exclusive of the 2.73 K microwave background radiation, as 

6xl0~ 15 _ 3 

Pu 7 er S cm I 1 ) 

ft 

where 0.05 < f t < 0.5 is the fraction that lies in the FIR/SMM range. The upper limit 
to ft comes from the UV/optical/near-infrared contribution; the lower limit comes from 
uncertainties, comprehensively reviewed by Dwek et al. (1998), and from the implications 
of the observed TeV gamma-ray flux from Mrk 501 (Coppi & Aharonian 1997; Stanev & 
Franceschini 1998). 



3. Contributions from Discrete Sources 

Using the Infrared Space Observatory, ISO, Taniguchi et al. (1997) found 15 sources 
at 5 to 8.5 /im, with a total flux of lmJy, in a 3' x 3' field, corresponding to a radiation 
density of p 7fim ~ 1.5 x 10~ 16 erg cm -3 . At 12 to 18/im, Altieri et al. (1998) derived a 
preliminary density of sources at fluxes above 50 /iJy at 2.5 x 10 -4 per square degree, 
yielding pi 5)im ^ 1.5 x 10~ 16 ergcm~ 3 as well. Kawara et al. (1998) conducted a deep survey 
over a 1.1 square degree field at 95 and 175 /im, respectively, and found 36 and 45 sources, 
all brighter than 150 mJy The corresponding radiation density in this ~ 140 /im-wide band 
at 120 /im is p 12 o^m ^ 1-5 x 10~ 16 erg cm" 3 . 

At 850 /im two deep surveys have been undertaken, both using the SCUBA array on 
the James Clerk Maxwell telescope in Hawaii. Hughes et al. (1998) reported five sources 
in a 5.6 square arc minute field for a total flux of 20mJy, or \I\ = 1.5 x 10~ 10 Wm 
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~ 2 sr _1 , while Barger et al. (1998) have found two sources with a combined flux of 8 mJy 
in two such fields of view. The respective energy densities at 850 /im are ~6x 1CT 17 and 
~ 1.3 x 10~ 17 ergcm~ 3 . 

4. Star Formation Rates, Metallicity, and Energy Production 

Madau et al. (1998) have made a persuasive case for estimating early star formation 
and energy generation rates, on the basis of luminosity of distant sources in the 0.15 to 
2.2 /im wavelength band. The top curve in Figure 1 shows their derived energy generation 
rates as a function of red shift. This approach is complementary to estimates of energy 
generation rates based on observed metallicities (Fall et al. 1996). 

Directly detected metallicity can underestimate the true heavy element and energy 
production. The fraction of processed material ejected from stars, f e j, and the energy per 
unit mass, e e j, generated in the process are only part of the story. We need to also consider 
the fraction of all baryonic mass, f(z), which has radiated away energy in producing heavy 
elements retained in the interior of stars. 

a. White dwarfs constitute a fraction f wd ~ 0.1 of the mass of all stars. Vassiliadis 
and Wood (1993, 1994) have traced the evolution of stars with initial masses ranging 
from 0.89 to 5 M Q , and find that the ejected mass contains hardly any elements heavier 
than helium, that the helium content of the ejected gas is only a few percent higher 
than the initial fraction, but that the final core masses of processed material range from 
~ 0.55 for the lower mass stars to ~ 0.9 M for stars at the high end of the range. We 
may, therefore, approximate the energy that was generated per unit white dwarf mass as a 
fraction e wd ~ 0.008 of the star's current mass. 

b. Neutron stars constitute a fraction f ns ~ 0.01 of all stellar mass, and have also 
converted a fraction e ~ 0.008 of this mass into energy. The neutron star's rotational 
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energy immediately after collapse into a body spinning with a millisecond period is also of 
order 0.2% of the mass-energy. This is eventually dissipated and radiated away. A typical 
neutron star with mass M ~ 1.4 M Q , thus will radiate away a fractional amount of energy 
per unit mass e ns ~ 0.01. 

c. The supernova explosion accompanying the formation of a neutron star provides the 
heavy elements and much of the non-primordial helium observed in stellar atmospheres and 
the interstellar medium. Chieffi et al. (1998) have recently compared estimates by various 
authors, of the total mass of heavy-element ejecta from supernovae, as contrasted to mass 
retained in the core. For a star whose initial mass was 25M Q the current consensus appears 
to favor ejection of a mass of approximately 4M in carbon, oxygen, neon, magnesium and 
silicon. In addition, the star also yields 9 M of 4 He, and a core mass of iron M Fe ~ 1.5 M Q . 
Consistent with the ratio of heavier elements to helium in supernova ejecta and stellar 
winds, we find that heavy elements constitute about 2% of the sun's mass, and that the 
sun's helium content is roughly 30% by mass, about 6% higher than the primordial helium 
abundance. The ratio of heavy elements ejected to those retained in neutron stars is also 
consistent with an overall neutron star complement amounting to ~ 1% of all stellar mass. 
We write f ej ~ 0.02, e e j ~ 0.01, respectively determined by solar heavy-element abundance 
and nuclear plus kinetic energies released in supernova explosions. 

d. Stellar mass black holes may radiate away an amount of energy comparable to 
neutron stars, but they appear to be rare, and probably do not contribute greatly to the 
overall electromagnetic energy constituting the extragalactic background. Neither the 
fractional mass in black holes at the centers of galaxies, nor the efficiency with which 
they have radiated away electromagnetic energy are currently well established, but their 
contribution to the metallicity and electromagnetic background appears to be modest. 

e. Finally, we need to also consider the class of low-mass stars that have never reached 
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the white dwarf stage, but have nevertheless contributed to the radiation density, though 
not to the observed metallicity. Their fraction by mass, fe m , depends on the stellar birth 
rate function. 

These different sources of energy suggest three separate estimates of contributions to the 
extragalactic background. The first is for a single burst, at red shift z, of rapidly evolving, 
massive, primordial, Population III objects with a fractional mass f e j(z) + f n8 {z). The 
second is for a continuous generation of massive stars that evolve primarily into neutron stars 
and white dwarfs. Here, the rate of formation at epoch z is f(z) = f e j(z) + fns(z) + f w d{z), 
which needs to be integrated over the appropriate red-shift interval Az. The last is for 
low-mass stars that may have formed at early times and have continued to shine at roughly 
constant luminosity for most of the history of the Universe. These cases are then confronted 
with current data. 



5. Must Most of the Energy Production have Occurred at Low Red Shifts zl 

The contribution to the radiation field by massive stars can be estimated from the 
dependence of the observed metallicity Z on red shift z. To lowest order, the data cited in 
Section 1 are consistent with a systematic metallicity decline by a factor of 10 for increasing 
red-shift intervals Az = 2 

Z/Z & ~ l(T z/2 ~ e~ L152 (2) 

5.1. A Single Star Burst at Red Shift z 

At epoch z, the baryonic mass density ps in a flat universe with deceleration parameter 
q Q = 0.5 is 

Mtfil + zf 

Pb Mb (3) 
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Here, is the baryon density parameter, H a is the Hubble constant today, and G is the 
gravitational constant. The fraction of mass that has been converted into heavy elements by 
epoch z, and the energy generated in this process is given by the stars that complete their 
evolution on time scales short compared to the age of the Universe and end up as supernova 
ejecta, neutron stars or black holes. We may approximate this by f(z) = al0~ z / 2 , where 
the value a ~ 0.1 is based on a solar abundance of heavy elements of ~ 2%, non-primordial 
helium of ~ 6%, and the complement of neutron stars that currently constitute roughly 1% 
of baryonic mass. Combining the nuclear and kinetic energies cited in section 4, we may 
take e ~ 0.01. 

If all Population III objects were formed in a single burst at epoch z, the metallicity of 
equation (2) and the corresponding radiation density p v observed today would be related by 

^,.(53^)^) a (4) 



or 



p v {z) = 1.5 x 10 15 f ^ rVf— ) (—) (—) (— W 2 ->/ 2 erg cm" 3 

i y! V50km s^^pc- 1 / vo.oiy \o.iJ vo.oiy \l + zj s 

(5) 

Comparison to equation (1) shows that a single burst, even as late as z = 2, fails to provide 
the requisite energy density for a value of Vt B as high as ~ 0.08. 



5.2. Continuous Formation of Massive Stars 

For continuous formation of massive stars from the earliest epochs, we note that the 
energy production rate is just the derivative of f(z) 

_ Vp. £ „ 1.15ea c e- 115z ~ 1.15ea c 10-*/ 2 (6) 

Here a c ~ 0.2, takes into account that roughly 10% of the baryonic mass is in white 
dwarfs, and assumes that the birth rate function ip does not greatly change over the aeons. 
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Integrating from z = 5 to the present, we have 

A,~6x ID- 4 ( ^ T ) Vik V M fA) erg cm -3 (7) 

P V50km s^Mpc" 1 / V0.01/ V0.27 V0.01/ 6 v ; 

This value is consistent with the observed background radiation even for Qb as low as 
~ 0.003, provided the FIR/SMM background is a fraction f t ^ 0.3 of the total. Since the 
metallicity rises sharply, all but ~ 4% of the light contributed to the background is of 
relatively recent origin, z <^ 2. 



5.3. Low-Mass Stars 

For low-mass stars, we simplify the calculation by assuming that they all have 
mass M < M , were formed at a single epoch, and have steadily radiated at constant 
luminosity, ever since. Taking a birth rate function ip oc M and a luminosity proportional 
to M 3 , we find the mean luminosity for stars of mass O.1M < M < M & to be 
(L/M) ~ O.O5L /M ~ 0.1 erg (g s) -1 . If a fraction f tm of all the baryons is in low mass 
stars, the generated radiation energy density increases at a rate p u (z), scaled down by the 
red shift. 



A^-Jt = { (L/M)(^S^-)f em (l + z)-Ht = -*to ( / 0Bfe ™ 
[1 + zY i W '\ 8ttG ) J K ' 2 °J (1 + z) 7 / 2 



where pob is today's baryon density and the age of the Universe at any given time is 
t ~ io(l + z)~ 3 / 2 . Stars that have been steadily shining since z max >> 1 contribute 

15 / (L/M) \/ H \ff tm \fn B \ _ 3 

to today's radiation density. This is significantly lower than the observed background. 
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5.4. Directly Observed Star Formation 

Finally, consider the observational star formation and energy generation rate given 
by Madau et al. (1998). The top curve of Figure 1 shows their data for the UV energy 
generation rate in a comoving volume of 1 Mpc 3 in units of erg Hz" 1 s" 1 on the left-hand 
scale. For the flat cosmological model with q = 0.5 that they assume, the figure also shows 
the corresponding energy generation rate per unit red-shift interval in units of erg Hz" 1 z -1 
shown on the right. Because time intervals At corresponding to unit red-shift intervals, 
Az = 1, diminish in proportion to (1 + z) 5 / 2 , this curve rapidly drops with increasing z. 
The red shift further lowers the contribution to the background by (1 + z), for a cumulative 
reduction by (1 + z) 7 ^ 2 . The integrated contribution to the background is ~ 40nWm" 2 sr" 1 
for an assumed UV/optical bandwidth Av ~ 3 x 10 15 Hz at the source. Recent suggestions 
that dust absorption in early galaxies would raise the data at high red shifts, to flatten the 
curve beyond the peak at z ~ 1.5 (Hughes et al. 1998; Pascarelle et al. 1998) hardly change 
estimates of the background radiation received today. Whatever the energy production rate 
at early times may be, the epochs at high red shifts do not last long, lead to an energy drain 
on the photons, and contribute little to the total energy density. However, for a c ~ 0.2, a 
high, massive-star-formation rate at early epochs would require an increase in Qb to keep 
the neucleosynthetic helium and heavier element abundances from exceeding their observed 
values. 

6. Discussion 

The FIR/SMM extragalactic radiation density is most readily understood in terms of 
energy generated by massive stars radiating primarily at recent epochs z < 2. Conversely, 
for currently estimated baryon densities, massive stars at high red shifts could not have 
generated the high levels of diffuse infrared background radiation observed with COBE. 
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Almost any star formation rate consistent with metallicities observed at different red shifts 
z leads to the same conclusion. Energy generation rates could have been high at early 
epochs z ^ 4, but the duration of such epochs is always brief except in inflationary Lemaitre 
universes, not considered here. 
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Foundation of Germany and the Max Planck Institute for Radioastronomy in Bonn. I 
thank both institutions for their warm hospitality. 
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Figure Caption 

Fig. 1 - Energy generation rates at different epochs, and their contribution to the present 
energy density. The curve with symbols is from Madau et al. (1998) and refers to the 
comoving UV energy generation rate, per second, for a flat universe with zero pressure 
and a deceleration parameter q Q = 0.5; here, the left-hand scale applies. The right-hand 
scale refers to the comoving energy generation rate per unit interval Az and to today's 
received radiation density. The energy generation and reception rates in both these sets 
of units remain invariant under red shifts, but the bandwidth diminishes as (1+z) leading 
to a reduced effective reception rate. These curves dramatically illustrate the negligible 
contribution that epochs at high red shifts make to the observed radiation density. To 
match the assumptions of Madau et al., all curves assume H Q = 50kms _1 Mpc -1 . 
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